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Abstract 

It is shown that in experiments on single molecule magnets (SMM's) in which transitions between 
two lowest spins states are induced by sweeping the applied magnetic field along the easy axis, the 
transitions are fully incoherent. Nuclear spins and the dipolar coupling of molecular spins are 
identified as the main sources of decoherence, and the form of the decoherence is calculated. The 
Landau-Zener-Stiickelberg (LZS) process is examined in light of this decoherence, and it is shown 
that the correct formula for the spin-flip probability is better given by a more recent formula of 
Kayanuma's than that of LZS. The two formulas are shown to be identical in the limit of rapid 
sweeps. An approximate way of incorporating the molecular spin dipole field into the rate equations 
for this process is developed. 
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I. INTRODUCTION 



A large number of molecular solids made from organic molecules containing magnetic 
ions have come to be known as single-molecule magnets (SMM's), and their magnetization 
dynamics has been studied intensively for over a decade now jj. The designation SMM 
comes about because the intermolecular magnetic interactions are much weaker than the 
intramolecular ones, yet one sees hystersis [2|, a phenomenon generally associated with 
ferromagnets in which the spins are strongly interacting. Of special interest, and the subject 
of this paper, is the study of low-temperature quantum tunneling between the two lowest 
Zeeman sublevels of one molecular spin (MS), since, then, processes such as phonon induced 
excitation or relaxation do not come into play jjj, and the dynamics is, a priori, purely 
quantum mechanical. 

The above conclusion is strongly reinforced by experiments in which the magnetization 



relaxes in the presence of a time- dependent magnetic fie 



where the two Zeeman levels are degenerate [J, 
classic Landau-Zener-Stiickelberg (LZS) process 



d which is swept through the value 
At first sight, this constitutes a 



9|, and the data appear to confirm this idea, 



especially in Fe 8 . The strongest check comes from the fact that the transition probability 
depends on the sweep rate over two and a half orders of magnitude in agreement with the 
LZS formula 5|. Further, the tunneling amplitude extracted by fitting to this formula agrees 
with direct numerical diagonalization of the single MS Hamiltonian. Most importantly, 
the matrix element so deduced varies with a static transverse magnetic field in oscillatory 
fashion 4], as required by the model Hamiltonian |l0| . 

It is, however, surprising that the LZS formula should be so well obeyed, since it is 
derived for isolated, noninteracting spins. The MS's in SMM's interact with many other 
degrees of freedom, and anything in the environment that can distiguish between the two 
tunneling states of the system will tend to suppress quantum tunneling and act as a source 
of decoherence. Phonons are an obvious such environment, but can be excluded by working 
at low enough temperatures. The remaining environment is that of the nuclear spins. These 



have been 
particles 



previously studied in connection with magnetization tunneling in small magnetic 

~\\ . In addition, one must also consider the other 



and in SMM's [jj, [JJ, 



MS's. The general picture that emerges from Refs. , ll4J is that nuclear spins give 

rise to incoherent transitions and the other MS's spin give rise to an additional magnetic 
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field that must be added on to the applied field in determining whether a given MS is at 
degeneracy or not. Other authors have adopted this point of view and studied these systems 
via Monte Carlo simulations [fj], Chapter 9 of Ref. [l| contains a good discussion of 
these and related points. 

Our purpose in this paper is to reexamine the decoherence from nuclear and molecular 
spins, especially the latter. In Fes, the dipole field from other MS's is about ten times 
larger than that due to the nuclear spins, so a priori, they should be a significant source 
of decoherence. It may at first sight be puzzling that the MS's which form the "system", 
can also behave as an "environment". The situation is analogous to how the electron- 
electron interaction in metals contributes to the electrical resistivity. In a model in which 
the MS's are coupled to each other, but not to any other degrees of freedom, the many-body 
(or many-spin) wave function of the MS's evolves coherently, yet the off-diagonal elements 
of the one-body (one-spin) density matrix can still decohere, i.e., decay with time. Since 
the magnetization is a sum of one-spin operators, such decay is relevant to its dynamics. 
Whether the model is adequate is a quantitative question depending on whether the omitted 
degrees of freedom are stronger or weaker decoherers than the ones considered. Thus, in 
metals at room temperature, phonon and impurity scattering are greater contributors to 
the resistivity than electron-electron scattering, and should not be omitted in a good model. 
The converse is true at very low temperatures in very pure samples (less than ~ 1 K in 
potassium, for example). 

The model we study is the following. Each magnetic molecule is taken to have a total 
spin S in its ground manifold, and to have two easy directions, ±z, separated by a barrier, 
Vb- It is assumed that other spin multiplets can be ignored at low temperatures, so that 
each molecule can be treated as a single spin of magntiude S. In zero external field, an 
isolated MS can tunnel between the m = ±S states. The corresponding energy splitting is 
denoted A. 

Next, the MS's are coupled to the nuclear spins (NS). Two broadly different types of 
couplings may be distinguished. If the magnetic ions have nuclei with nonzero magnetic 
moments, the contact hyperfine interaction between an ion and its own nucleus must be 
considered. The corresponding energy scale is 1-10 mK. The second is the dipolar coupling 
between the MS's and other nuclear spins, with an energy scale E^n ~ 1 mK for close by 
nuclei. (The suffixes 'd' and 'n' stand for 'dipole' and 'nuclear', respectively.) We shall 
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assume that E dn ^> A, as is the case in Feg. 

In addition, different MS's are coupled via the dipole-dipole interaction, which is taken 
to have a scale ~ E^ m for nearest neighbours. (The suffixes 'd' and 'm' stand for 'dipole' 
and 'molecular', respectively.) There is clear separation of energy scales: Vb S> E dm ^> A. 
This is a good description of many SMM's. In Fe 8 , e.g., Vb ~ 20 K, E dm ~ 0.1 K, and 
A ~ 10~ 7 -10~ 8 K. Stray and dipolar magnetic fields along x and y are unimportant since 
they are not large enough to give any significant mixing of the m = ±S states with the 
higher Zeeman states, and they affect A only weakly. Along z on the other hand, such 
fields are very important, since they move MS's off resonance. Under these conditions, each 
MS may be replaced by a pseudospin with spin-1/2 with the | |, j) states representing the 
m = ±S states of the true spin. 

The plan of the paper is as follows. We calculate the decoherence from nuclear and 
molecular spins in Sees. [IT] and 1111} respectively, pushing various details of the calculations 
to the Appendices. In Sec. IIVI we consider the two environments together. In Sec. |V] we 
consider the implications of the decoherence for the LZS process. We find that although 
the tunneling is indeed incoherent, the net spin-flip probability in a single LZS sweep is 
remarkably insensitive to the details of the decoherence mechanism. In a simple model 
where the dipole field due to the other MS's is omitted, the probability turns out to be 
given exactly b y K ayanuma's formula for a spin coupled to an oscillator bath in the strong 
damping limit 13] • In the limit of high field sweep rate this formula agrees precisely with 
the LZS formula. This explains why the experiments appear to be in accord with the 
LZS scenario. We also consider a better approximation where the dipole field is included 
in a macroscopically averaged way. This approximation improves the agreement with the 
experiments by Wernsdorfer et al. 0, 6]. 



II. MODEL FOR NUCLEAR SPIN ENVIRONMENT 

As our first model, we consider a single molecular spin interacting with the nuclear spins 
via the dipolar coupling. Hyperfine and transferred hyperfine interactions are not explicitly 
included, although in the end they are unlikely to have qualitatively different effects, and 
only to lead to a modification of the parameter W introduced below. We assume that all 
nuclear spins have spin 1/2, and neglect the local magnetic field i?i oc at the nuclear site. 
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This is a good assumption if H\ oc <C ksT/fin, where /i n is the nuclear magnetic moment. 
This is indeed so since fc^T '/ '/x n ~ 10 T at 10 mK. The dipolar coupling between nuclear 
spins can be neglected for the same reason. With these assumptions, our Hamiltonian is 

H mn = ]r(Acr 0x + ea 0z ) + Yj l a oz(?iz - 3er 02 cos W 5< • r 0i ]. (2.1) 

1 i r 0i 

Here, % labels the different nuclear spins, c?q and a, denote the Pauli spin matrices for the 
MS and the ith nuclear spin, r j is the position of the ith NS relative to the MS, r 0i = \r 0i \, 
and cos^oi = z ■ roi/roj. Further, a is the characteristic distance from the MS to the nearest 
NS. We expect a ~ 1-2 A for any SMM. Finally, we have included an energy bias e between 
the | |) and | |) states of the MS, which could arise from an external magnetic field. The 
suffixes in H mn stand for 'molecular' and 'nuclear'. 

We now suppose that at time t = the MS is in the state 1 1), and that every NS is in a 
completely disordered state described by the density matrix 1/2. Again, this assumption is 
well justified at the temperatures at which experiments have been carried out so far. The 
quantity of interest is the probability, P(t), that the MS will be in the state | j) irrespective 
of the NS state. 

Even for this simple model, an exact calculation of P(t) is not possible (but see below). 
We therefore turn to the approximate methods described in Sec. Ill A-D of Ref. 18]. 
We cannot assume that the damping is weak, or that the NS's are fast compared to the 
MS's. A "golden rule" approach is still fruitful, however, as A is the smallest energy scale 
in the problem. Moreover, the validity of this approach can be self-consistently checked. 
Second-order perturbation theory yields 

P(t) = ^ fdh fdhe^-^HF^h), (2.2) 
4 jo Jo 



where 

FJti.tn) = ^ 



Fi(ti,t 2 ) = -Tii [e iWl+tl e m '- {tl - t2) e - im+t2 } , (2.3) 



with 

H i± = ± Edn 3 a [<j iz -3 cos 9 0i ai ■ r Qi \. (2.4) 

The quantity Fi is the contribution of the zth environmental spin to Feynman's influence 
functional evaluated for a particular pair of forward and backward paths of the "system" 
spin, namely, that in which this spin flips from up to down at time t\ on the forward path, 
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and time t% on the backward path. We therefore refer to Fi as the (environmental) influence 
factor or function. 

The trace in Eq. (12. 3p is easy to evaluate. Defining 

ti2 = h-t 2 . (2.5) 

x oi — r oi ' x > etc., and the vector 

h, = Edr f (-3z 0l x ou -3z 0i y 0i , r 2 m - 3a&), (2.6) 

we have 

F i (t 1 ,t 2 ) = cos2/i i ti 2 . (2.7) 

Now, 

^ = ^!(l + 3cos 2 ^) 1/2 , (2.8) 

so hi ~ Edn for the nearest NS, and drops as 1/r 3 for more distant ones. Thus, for £i 2 > Ej^ 
the different F^s have random signs, and since they can not exceed 1 in magnitude, they 
essentially multiply out to zero. We conclude that phase coherence is lost on the time scale 
t c ~ EJn, and for t t c , we get incoherent tunneling. For such times, we can approximate 

n^(ti,t2)^exp(-2^/ ij 2 4). (2.9) 

i i 

Further, in the double integral in Eq. (12.21) . we may introduce sum and difference variables 
i = {fi +t 2 )/2 and r = t 12 . The integral over r is essentially independent of i, and its limits 
may be extended to ±oo. The i integral then gives an overall factor of t, yielding 

P(t) ^ T n t, (2.10) 



where, with, 



W 2 = 4j2hl (2.11) 



i 



1 . o r°° , 1 H/U V27T A 2 



-A 2 dre^e-^ = l^^ e -?l™\ ( 2 .12) 

We may estimate W by replacing the sum in Eq. (12.111) by an integral, taking a uniform 
density of nuclear spins equal to 1/a 3 outside a sphere of radius a. Since 

h 2 = ^-(l + 3cos 2 e 0l ), (2.13) 



r, 



6 

6 



W 2 ^4Ela e l ^(1+3 cos 2 61)1 (2.14) 



> a a 3 r 6 
= "fEl. (2.15) 

In fact, the integral estimates the contribution of the nearest neighbors rather poorly, and for 
the simple, body-centered, and face-centered cubic lattices, the number multiplying E\ n is 



67.2, 98.0, and 116, respectively [191 ] . Thus, in order of magnitude, we may take W ~ lOE^n 
for any magnetic molecular solid. It should be noted that for a fixed bias e, the rate T n goes 
up with increasing Ed n , as long as e > W. The converse is true for the very small number 
of MS's on which the bias is small, e < W. 

The result ( 12. 10[) is essentially a Fermi golden rule rate, and is limited to t -C T" 1 . For 
longer times, a formal answer can be obtained as follows [201 ] . We can write 

H mn = ^A-a , (2.16) 

where 

A = Ax + (e + 2^hj • <7;)z. (2.17) 

i 

Thus, A is an operator with respect to the bath spins. With the understanding that these 
must be traced over, we get 

iA 

(i \ e - lHmnt \ t) = -— sin \kt. (2.18) 

Thus, 

^) = A 2 n^(lsin 2 l), (2.19) 

where tr^ indicates a trace over the ith NS. To perform this trace we take the quantization 
axis for it to be parallel to hj. This means that the variable 

B n = 2j2hiSi (2.20) 



takes on all possible values obtained by letting each be +1 or —1 independently [2l|]. That 
is to say, B n is a stochastic variable with some probability distribution, P(B n ), and the spin 
flip probability is obtained by averaging over this distribution: 

/oo A 2 
^ — sm\\M)P{B n )dB n , (2.21) 
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with 



1/2 



A= [A 2 + (e + B n ) 2 y". (2.22) 

To proceed further, we need the form of P(B n ). We find this approximately by arguing 
that because of the law of large numbers B n is a Gaussian with a variance W 2 , i.e., 

1 \ V2 



P(B ri 



-BI/2W 2 



(2.23) 



,2ttW 2 . 

We can do somewhat better by looking at the moments of B n . We clearly have (B 2 ) = W 2 , 
but 



(B*) = 3(B 2 n ) z -32j:ht 



(2.24) 



Thus the fourth moment is less than what it is for a Gaussian (negative kurtosis), and the 
distribution has less weight in the wings than a Gaussian. We shall see that that the detailed 
form of P(B n ) is not too important, and for our purposes, Eq. (12.231) is good enough. 

For W~ l < t < A -1 , we may evaluate P(t) by replacing A by (B n +e). (This replacement 
is no longer valid when At > 1, for then the phase of sin 2 (At/2) is significantly altered by 
throwing away A.) Then by the usual textbook argument for Fermi's golden rule, 

sin 2 ((S„ + e)t/2) 2vrt 



(B n + e) 2 

The integral for P(t) is then trivial, and yields 

-/2^A 2 



P(t) 



—5(B n + e). 



-e 2 /2W< 



t, 



-e, and we may put B r , 



4 W 

which is the same as before. 

For At > 1, the integral is dominated by B n ', 
Gaussian factor, yielding 

A 2 _ t 2 /2W 2 r°° 1 - cos(VA 2 + b 2 t) 
J-oo A 2 + b 2 



(2.25) 



(2.26) 



-e in the 



P(t) 



db 



J (z)dz 



(2.27) 



8W v J At 

where b = B n + e, and we used Ref. {22] in the last step. Using the asymptotic behaviour of 
the Bessel function, we find that for At ^> 1, 

fWA 



P(t) 



-e 2 /2W 2 



8 W 



1 - W^-sinfAt - - 
V vrAt V 4 



^2.28) 



The important point is that even for e = 0, the nuclear spin environment impedes the spin 
from flipping appreciably, and the net flip probability is only of order A/W. 
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III. MODEL FOR MOLECULAR SPIN ENVIRONMENT 

For our second model, we consider only the dipolar coupling between MS's, and ignore 
the nuclear spins altogether. Let us denote the energy scale of the mutual dipole-dipole 
interaction between MS's by E^m [see Eq. ( 13 .3p below for the exact definition]. Since E^m 3> 
Edni we may a priori expect decoherence by the mutual interaction to be much greater than 
that by the interaction with NS's. This model is studied in an attempt to investigate this 
point. 

In terms of the Pauli matrices, the Hamiltonian for interacting MS's can be written as 

= qXX^ 7 ** + e i a iz) + ■Z^KijVizVjz- (3.1) 

1 i 1 i<j 

Here i and j label the different spins, x and z denote the axes, 6{ is the bias field on spin i 
that moves it off-resonance, and is the dipolar coupling. 

Let us now focus on one MS, which we shall call the system, and label it with a suffix 
0. This is prepared in the | |) state at time t = 0, and the other spins, which we call 
the bath, are prepared in a density matrix pb- Let P(t) denote the probability that the 
system spin is in the state | j.) at a later time t irrespective of the state of the bath. For 
an isolated spin, P(t) = sin 2 (At/2). If decoherence is weak, we expect the oscillations to 
be weakly damped, and if it is strong, we expect a decay without any oscillation. Indeed, 
these qualitative behaviours define what we mean by weak and strong decoherence. Since 
the dipole interaction is long-ranged, we anticipate that the decoherence might depend on 
the spatial position of spin in the sample, especially if pe corresponds to a fully or nearly 
fully polarized bath, but otherwise there is nothing special about its choice. 

The calculation of P(t) for the model (13. ip appears daunting because of the couplings 
between the bath spins. We therefore consider a simpler model 

7-tmm = -(Aooz + ecroz) + 7: ^2(Aa ix + iz ) + - Kia 0z a iz . (3.2) 

(Both suffixes in 7i mm stand for 'molecular'.) The dipolar couplings between the bath 
spins are now replaced by a distribution of dipole fields by treating the bias energies q as 
independent random variables, distributed on the scale Ed m . The calculation of P(t) should 
include an ensemble average over this distribution. The coupling Ki between spin and 
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spin i of the bath is, however, retained as is, and is, explicitly, 

2E dm a 3 



Ki 



1 - 3 cos e 2 0i ) 



(3.3) 



0; 



Here, a is the nearest neighbour distance, is the distance from spin to spin i, and #oi 
is the angle the line joining them makes with the z axis. Finally, e is an additional bias on 
spin 0, due to an external field, for example. 

For purposes of explicit calculation, we shall take the probability density of the biases 
to be Gaussian, 



i 



2 / 2E i 



(3.4) 



where Ei, ~ Ed m . Dipolar field distributions in Fe§ have been measured by Ohm, Sangre- 
gorio, and Paulsen [23], and by Wernsdorfer et al. 24]. They have also been inferred from 
linewidth measurements in optical spectroscopy by Mukhin et al. [251 ] . The assumption 
of a Gaussian form is consistent with these measurements. Berkov has given theoretical 
and Monte Carlo arguments for a Gaussian distribution in a system of dense interacting 



dipoles 



26] . We shall see, nevertheless, that the detailed form of this distribution is not 



physically important for us. 

Even the model ( 13. 21) cannot be treated exactly. It is again seen that the weak coupling 
approximation is totally invalid, and adiabatic renormalization is inapplicable since the 
bath and system spins move on the same time scale. The golden rule is still good, however. 
Second-order perturbation theory in A yields 

.2 



A 2 r l /•* , 
Pit) = — dt x dt 2 e te{tl - t2) F, 
4 Jo Jo 



where 



with 



Tr, 



n 



PbU 



JHi+h -iHi-(ti-t 2 ) -iH i+ t 2 



i± 



I (A 

9 V 



a ix + (ei ± Ki)o L 



(3.5) 



(3.6) 



(3.7) 



The choice of pe demands some care. It would now be incorrect to take pe — 2~ N,n , where 
N m is the number of MS's since these spins do not equilibrate between the | f) an d | |) 
states on a time scale short compared to A -1 . Instead we choose each spin to be in a definite 
state, either | f) or | j). ( m the language of statistical mechanics, the bath is in a state 
of quenched disorder.) This then means that in principle we have to calculate the second 
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order influence function for every configuration of MS's separately. In practice, this is not 
so, and we shall see that the functions for up and down spins differ only by phases. When 
these phases are added together for all the MS's in the bath, they will reproduce exactly the 
effect of the local dipole field at spin 0. This field is dependent on the MS configuration, 
but except for special configurations such as all or nearly all MS's polarized in the same 
direction, we can treat it statistically as a field with an rms value of order Edm- 

In equations, the above means that if we specify the spin configuration by giving Sj = 
(a iz ) = ±1, then 

PB = Y[pi] Pi = 7^(1 + s i°iz)- (3.8) 

i ~ 

Accordingly, F factorizes into a product of factors, one for each bath MS. If the zth spin is 
"up" , this factor is 

F { = (t | e ^+*i e - i ^-(*i-fe) e -^+*2| (3 9) 

If the spin is "down", Fi is given by the expectation value of the same operator in the | [) 
state. The calculation of these influence factors is lengthy, and is presented in Appendix lAl 
We find that 

F i ^e iSiK ^{l-r ]i ), (3.10) 
with rji given by Eq. (1A46I) with the addition of a suffix i to K and Q±, i = (ti + 12)/2, and 

tl2 = ti — ^2- 

We call the quantity rji the mismatch, since it arises from a difference in the time evolution 
of the ith environmental spin in response to different paths taken by the system spin. The 
derivation in Appendix IA1 shows that < < 1, vanishing only when t 12 = [27]. Hence 
we may put 1 — rji « e~ Vi , leading to 

P(t) = ^ f dt x f dt 2 e ieTtl2 e-^i r, % (3.11) 
4 jo Jo 

where 

e T = e + Y,K l s i . (3.12) 

i 

This is the total bias that the spin at sees including the dipole field of the other MS's. Its 
value is of order Edm except for special spin configurations. 

We show in Appendix [B] that for |ti 2 | > E^, and A" 1 < t < EjJA, 

J2vi-lmA\ti 2 \, (3.13) 
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where 7 m is a constant of order unity. We have also evaluated this sum numerically, as 
described in Appendix IB 21 This work shows that the form (13. 13[) is good even for A|£i2| ~ 1. 
Employing it in Eq. (13. lip , we get 



A 2 

P(t) = —Re 



I l — e -(7mA-ie T )i 



(3.14) 



_ (7 m A - ie T ) (7mA - ie T ) 2 
Thus, Pit) displays damped oscillations about a slowly rising mean. The time scale of the 
decoherence is A -1 , which is comparable to the time scale of the oscillations when the total 
bias, €t, is zero. The amplitude of the oscillations is ~ A 2 /e|, if the bias is large. For 
t A -1 , we obtain 

p(t)~r m t, (3.15) 

with 

1 7 A 3 

r ™ = 5^#T4' (3 ' 16) 

This quantity may be interpreted as an average rate at which the spin flips. If the net bias 
is large (s> A), this rate is 7 m A 3 /2e 2 -., while if the bias is zero, it is much larger, A/27 m . 
(The amplitude of the oscillations is also very small when the bias is large.) It is interesting 
that the zero-bias rate is proportional to A and not to A 2 as might be expected from a naive 
application of the golden rule; this is because the decoherence time scale is also set by A. 



IV. COMBINED NUCLEAR AND MOLECULAR SPIN ENVIRONMENTS 

Let us now consider both environments together. The combined influence factor is the 
product of the influence factors for each separate environment, leading to 

P(t) = — t dh /' dt 2 e i£T * 12 e- 7mA| * 12| e- w2 *"/ 2 . (4.1) 
4 jo Jo 

If, as is generally the case, W ~ E dn ^> A, the integrals may be evaluated as in Sec. [Ill We 
once again get P(t) ~ Ft, with 

A 2 poo 

T = ^J dte^e-^e-^ 2 ' 2 . (4.2) 
In general this integral leads to an error function, but if E dn 3> A, it simplifies, and we get 

r = ^£-*""- < 4 - 3 > 

This is of the same form as T n , and the main effect of the molecular spins is to change the 
bias field. 
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V. QUASISTATIC MODEL OF FIELD SWEEPS 



We have seen in the previous sections that the molecular spin relaxes incoherently from 
| |) to | |). There may in addition be some vestige of the coherent oscillations, but these 
decay because of the coupling to nuclear spins and to other molecular spins. The decay 
time scales due to these two couplings are Ej^ and A -1 respectively, and the former is 
the relevant one since it is so much shorter. If the externally applied field is swept slowly 
enough that the bias on any one spin changes by much less than Ea n in a time E^, that 
is, if €t -C E% n , then it is a good approximation to neglect the off-diagonal elements of the 
density matrix, and to write simple rate equations for the diagonal elements. If we denote 
the probability for a particular molecular spin to be in the | f ) or | |) states by pj and 
we have, 

^ = r(e r (t))(a - pt) = r(6r(f))(l - 2p T ). (5.1) 

where the rate V has been allowed to vary with time through its dependence on the bias. 
Let the spin state be | j) at t = — oo. Then, Eq. (15. ip is easily integrated to yield 

1 



exp (-2 J' T(e T (t'))dt'^ 



In particular, the probability for the spin to flip is given by 

, , 1 



Pf =p^oo) = - 



1 - exp (-2 r Y{e T {t'))dt' 

\ J — CO 



(5.2) 



(5.3) 



It is interesting to analyze the spin flip probabilty neglecting the contribution of the other 
molecular spins to the bias. That is, we take er{t) to be e a (t), the applied bias field. Further, 
as in the standard LZS protocol, we take e a to be a constant. Such an analysis would be 
directly applicable to a situation in which the molecular spins were very dilute and E^m was 
smaller than Ed n - Since we chose oo) = 0, we must take the bias field to be swept from 
large positive to large negative values and the integral in Eq. (15. 3p becomes 

/oo I POO 

T(e a (t))dt = — T(e a )de a . (5.4) 
-oo £ a J -co 

Equation (14.21) now yields (writing e a for et) 

^2 „oo 



/oo /\* rco rco „ „ 

T(e a )de a = dej dt e - V^'e"^ 2 / 2 

-oo 4 J —oo J oo 

A 2 poo 

^ J dte-^e- wH2 ' 2 2it8{t) 



A 

ttA 2 



(5.5) 
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Hence, 



p / = i(l-e- A2 /H (5-6) 



This is the same as Kayanuma's result 17| for a spin coupled to an oscillator bath in the 
strong damping limit. Our derivation shows that this result is valid more generally whenever 
the transitions are so incoherent as to allow for rate equations. The striking fact that the 
details of the decoherence mechanism drop out of the final result can also be understood. 
If the decoherence is large, the rate T is small, but the spin can flip over a larger energy 
interval around the crossing, i.e., over a larger range of bias energy. For pure nuclear spin 
decoherence, T ~ A 2 /W, but the crossing region is broadened to a width ~ W. For pure 
molecular spin decoherence, T ~ A, and the crossing region is also of width ~ A. 
The quasistatic result (15. 6ft should be compared with the LZS spin-flip probability, 

P/,Lzs = (l-e^ A2/2|ea| ). (5.7) 

In the fast-sweep limit, i.e., with |e a | ^> A 2 , pf <C 1, and the two results are identical, 

ttA 2 , . 

Pf =P/,lzs - rp-j-- (5.8) 



This remarkable resu 



t has very interesting implications for the experiments by Wernsdorfer 
(|. It has always been a surprise that the data in these experiments 



and colleagues fly, 

agree with the LZS formula, even in the fast sweep limit. After all, the LZS formula is 
derived for a single noninteracting spin, and the spins in Feg are not noninteracting, and are 
subject to strong and rapidly fluctuating fields from the NS's and possibly the MS's. Indeed, 
it is the systematics of the agreement with the LZS formula that has been used to argue that 
one can extract the underlying tunneling matrix element from the incoherent relaxation of 
the net magnetization in a swept external field. Equation (15. 8p provides an explanation of 



this fact. It also means, in a stroke of luck, that the analysis of Ref. 28J continues to be 
valid. 

In the slow-sweep limit, on the other hand, P/,lzs — 1) while pf ~ 1/2. This means that 
if we continue to infer a tunneling matrix element, A in f, by fitting the flip probability to an 
LZS form, we have 

Af nf (e) = -^m [t±±- . (5.9) 

We plot this in Fig. [TJ which should be compared with Fig. 7 of Ref. |5J . Although our plot 
is qualitatively similar, it does not agree in detail. In particular, the experimentally inferred 
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splitting drops more rapidly with e once e < 50A 2 than our model shows. Neverthless, 
the general trend indicates that we have captured some of the essential physics. On the 
other hand, this simple formula does not contain the experimentally seen dependence of the 
splitting as inferred from the fast-sweep data on the nuclear spin coupling [6]. 

To prevent misunderstanding, we note here that we have described the sweep as fast or 
slow depending on the ratio |e a |/A 2 . However, because A <C Ed n , even if |e a | A 2 , it is 
possible to satisfy |e | -C Ej n , the condition for the quasistatic treatment to apply. 

Let us now ask how to include the effect of the mutual dipole field as the external magnetic 
field is swept. The picture that emerges is that each MS flips at a rate that depends on the 
bias field seen by it. If at time t, the MS configuration is set of Ising spin variables {s^}, 
and the net bias on the ith spin is e^, then at a short time At later, 

Si — > —Si with probability T(eiT(t))At. (5.10) 

The spins and the dipole fields then become a complicated coupled stochastic process. As 



noted in Ref. [1], this is a Glauber process with the difference that the flipping rate depends 
on the long ranged dipole field. Monte Carlo studies of such processes have been performed 

n n 

by Cuccoli et al. [15J, and by Fernandez and Alonso [16(. It would be interesting to conduct 
similar studies in a swept field with the rates found by us. 

Here, we consider a simpler way to incorporate the dipole field in the rate equation 
(Eq. (15. ip ) in an average way that ignores its site to site variation, through the macroscopic 
demagnetization field. To forestall confusion, it pays to recall the distinction between B, H, 
and the contribution of the demagnetization field to the latter. We work in the Gaussian 
system of units. Let H a be the applied magnetic field, i.e., the field that a solenoid wound 
around the sample would produce if the sample were not there. Let M be the magnetization, 
i.e., the magnetic dipole moment density, and let Hdemag be the demagnetization field, i.e., 
the field produced by a volume charge density V • M and a surface charge density M ■ n, 
where n is the outward normal at the surface of the sample. The field to which an individual 
MS responds is the induction 

B = H a + 4vrM + H demag , (5.11) 
through a term in the Hamiltonian, 

Hbias = -gfi B S op ■ B. (5.12) 
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Here, S op is the operator for the total spin of the molecule in question. Since, as we have 
argued, the MS's behave as essentially classical variables with only a z component, the total 
bias on this MS is 

e T = -2 9f M B SB z . (5.13) 

Henceforth we will write just H a and M for H a z and M z . For simplicity we will ignore 
the spatial inhomogeneity of M and Hdemag, as well as the tensorial character of their 
proportionality, and write 

4ttM + H demag = oM. (5.14) 

The constant a is shape dependent: it would be 8n/3 for a perfectly uniformly magnetized 
sphere, 4tt for a thin long rod parallel to H a , and for a thin flat disc normal to H a . 
With the above definitions, the bias is given by 

€ T = -2gfi B S(H a - angfi B S(l - 2p T )), (5.15) 

where n is the number density of MS's. Hence, n(gfi B S) 2 ~ E dm . Let us again take 
p^(-oo) = and H a < at t — > — oo, so i a < 0. Adjusting the zero of time, and absorbing 
another constant of order unity in a, we get 

e T (t) =e a t-4aE dmP i(t), (5.16) 

and feed this into the rate equation (15.11) . The resulting differential equation for p^ is 
rfp T _ v 7 ^ A 2 



1 - 2p T )exp 



1 / 4aE dm 

-2 r + ^r v ^ 



(5.17) 



du 4 |e a | 

where u = \e a \t/W. This equation can also be formally solved as before, by treating r(er(t)) 
as a known function of time. The solution is then again given by Eq. (I5.2p . but since er(i) 
depends on p^(t), it is now in the form of an integral equation. We have found it simpler 
to integrate the differential equation numerically for different values of e. The results are 
shown in Fig. [TJ The qualitative agreement with Ref. Q| is improved, although we cannot 
make a direct comparison because of uncertainty in the ratio aE dm /W. 

Finally, let us return to the point that for slow sweeps, Pt(oo) ~ 1 for coherent LZS 
sweeps, while p-f(oo) ~ | for incoherent sweeps, both from Kayanuma's formula (15. 6p or the 



formula (15.171) which includes MS dipolar fields 29|. This means that starting from a sample 



with a saturated magnetization —Mo, we are arguing that the final magnetization will be 
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and not Mp if the field is swept slowly. And indeed, studies of the Mn 4 |30j and Mn i2 wheel 



SMM's [3l| show just such behaviour. In Ref. |30j, it is found that the final magnetization 
is zero for slow sweeps (see Fig. 8a, 8b, and 8c there), and is fully reversed only for ultra 
slow sweeps (Fig. 8d and 8h). For sweep rates in between, and for inverse LZS sweeps (Fig. 
8c, 8f, and 8g), the final magnetization is not zero, but is not completely reversed either. (See 



also Fig. la of Ref. [3l|. This supports the conclusion that the transitions are incoherent. 
It also means that to fully explain the ultra-slow sweep and the inverse LZS sweep data, 
one must have a mechanism for the spin to relax from the higher energy state to the lower 
energy one even (but not vice versa) when the bias is much larger than Ed n - One possibility 
is to have a second order Fermi golden rule process in which (assuming the spin is 10), the 
spin tunnels from the m = — 10 to an m = 9 or m = 8 virtual state followed by a transition 
to the m = 10 state with the emission of a phonon. We shall address this issue further in a 
separate publication. 

VI. DISCUSSION 

We have considered the transitions in a swept field in the presence of nuclear and molec- 
ular spin decoherence. Our qualitative conclusions regarding the former are in accord with 
those of Refs. Ij], Q, but the quantitative form of the decoherence is different. Sim- 
ilarly, with regard to the molecular spins, we agree with them and other authors 
that their main effect is to add an essentially c-number contribution to the bias field on any 
given MS. However, we believe that this conclusion was not foregone, and that our treatment 
gives a proper justification for neglecting the additional decoherent effect of these degrees 
of freedom. 

The quasistatic approximation enables us to answer the question posed at the start, viz., 
why the LZS formula appears to describe the swept field experiments so well. We find 
that this is not because the transitions are coherent, but because the effective width of the 
crossing and the incoherent spin-flip rate vary inversely, leading to a fortuitous cancellation. 
It remains an open question to study the stochastic variation of the bias field, and thus 
understand this process even better. 
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APPENDIX A: CALCULATION OF SINGLE-SPIN INFLUENCE FACTORS 
FOR MOLECULAR SPINS 

In this Appendix, we calculate the factor Fi in Eq. (13.91) . To save writing, we omit the 
index i henceforth. Let us denote the influence factor by or F^ when the environmenal 
spin is up or down respectively. We further abbreviate 

e± = e ± K, (Al) 
fi| = A 2 + e|, (A2) 

and £12 = ti — t 2 as before. 

Let us find F^ first. As a first approximation, we argue that because A is much smaller 
than the typical value of e or K, we may neglect it altogether. This yields 

F^e lKtl2 . (A3) 

This approximation is too crude. It implies \F^\ = 1, which is the maximum possible value 
it can have. Decoherence arises precisely from the fact that \F^\ < 1 because [H+,T~C-] 7^ 0. 
It is important to find the departure from unity. With this in mind, let us write 

F t = c*(l-f7), (A4) 

where is a real phase, and 77 is another real quantity that we have refered to as the 
mismatch. Our crude calculation shows that 77 is small and ~ Kt\2- 

Before calculating r\ more carefully, let us relate to Fj. Since | |) = —icr y \ f), we can 
write 

Fi = (T W y e H ^a y a y e- m -^-^a y a y e- m ^a y \ |), (A5) 
and since a y anticommutes with 7i±, this can be transformed to 

F l = (T \ e -i-H + t le iH-(ti-t2) e m + t 2 \ T ^ (A6) 
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which is the same expression as i*j with the signs of t\ and t 2 reversed. That is, 

F l (t 1 ,t 2 )=F 1 (-t 1 ,-t 2 ). (A7) 
In particular, we shall see that the mismatch for F± is the same as that for F^, so that 

F^t^^e-^il-rj). (A8) 
To find the mismatch more accurately, we write 

F T = ( ni |n 2 ), (A9) 

where 

|n a ) = e in-t ae -m + t a \ ^ a = l,2. (A10) 

The notation in this equation exploits the fact that every pure state of a spin- 1/2 system can 
be written as a spin-coherent-state, i.e. a state with maximal spin projection along some 
direction in space. Thus, the states defined in Eq. ( lAlOj) have maximal spin projections 
along directions fii and 112. These directions remain to be found. Of course, the states also 
have phases which also need to be found. 

Let us now view Eq. flAlOl) in terms of two rotations applied to the state | f). Since 
A <^ e, these rotations are both about directions very close to z. Accordingly, fii and n2 
are also very close to z, and we may write 

n az ^ 1 - \n 2 a± , (All) 

where n a ± is the component of n a perpendicular to z. Now, since 

\(ni\n 2 )\ 2 = -(l + n 1 -n 2 ) (A12) 



for spin-1/2 coherent states, we may write 

I^Tl 2 « \ (2 " \n\ ± - \n 2 21 _ + n 1± ■ n 2± ) (A13) 
= l-i(n 1± -n 2± ) 2 . (A14) 

Taking the square root, and recalling the definition of the mismatch, we get 

V ~ kni-L - n 2± ) 2 . (A15) 

o 
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The problem is thus to find n a ±. We have not been able to find any simple way to 
do this except by explicit expansion and multiplication of the exponentiated operators in 
Eq. OAlOp . The resulting trigonometric expressions can be made somewhat easier to handle 
if we introduce the abbreviations 

1± = Otti/2, 6 2± = Ofcta/2. (A16) 

With these, we may write 

-lH+t 2 



where 



Similarly, 



with 



Then, 



with 



1 

c = cos # 2+ , c = - — — sin6 l 2 + (A, 0, e+). (A18) 

i Z_i_ 



3 *w-*a = b + b • a, (A19) 



b o = cos0 2 -, b = ^-sinfl 2 -(A,0,e_). (A20) 



|n 2 ) = {b + b ■ a)(c + c • <?)| T) = A 2T | T) + A n \ |), (A21) 



^2 T = {bo + b z )(c + c z ) + b x c x , (A22) 
A 2 .i = (bo-b z )c x + b x {co + c z ). (A23) 

In terms of these quantities, we have 

n 2 + = n 2x + in 2y 
= (n 2 |a + |n 2 ) 

= 2A^A U . (A24) 

By comparing and real and imaginary parts of both sides, we obtain n 2x and n 2y . We now 
note that the quantities A%\ and A 2 i consist of various terms oscillating at the sums and 
differences of the frequencies Q±/2. Since A Q± for all but very distant (and therefore 
very weakly coupled) MS's , we may expand the amplitudes of these oscillatory factors in 
powers of A. Using the results 



e + A^ , k N 

— ~ 1 " W {A25) 
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etc., we obtain 



c + c z = e- lf?2 + + 0(A) 2 , 
b ±b z = e ±l9 *- + 0(A) 2 . 



(A26) 
(A27) 



Therefore, 



A 2 -\ = e^ 2 — e2+) + 0(Af 



M = 



and 



A - a A fl 

-i— sm9 2+ e- td2 - + i— sin0 2 _e~ if?2 +, 



(A28) 
(A29) 



(A30) 



From this expression, we can get n 2:r and n 2y by taking and real and imaginary parts. 



n 2l 



-2 — sin 2+ sin(20 2 _ - 2+ ) + 2— sin 2 2 _ 



A 

A 2AK 



= - — [cos2(0 2+ - 02-) - cos20 2 _] + A[i _ C os20 2 _ 

cos2^ 2 _ - -A- cos2 (^ 2+ -e 2 -), 



where in the last line we have used the result 

1 1 



IK 



In the same way, we have 



(A31) 
(A32) 
(A33) 

(A34) 



n 2y = -2^- sin 6 2+ cos(20 2 _ - 02+) + 2— sin 2 _ cos 2 _ 



--A-[ sin2 (# 2+ - 2 -) + sin 20 2 -] + A s i n 202_ 
i i I— 

2AK A 

sin20 2 _-— sin 2 (0 2+ - 2 -). 



(A35) 
(A36) 
(A37) 



For n lx and n ly , we simply change the suffix 2 in 2 ± from 2 to 1. We then have 

nix ~n 2x = AA^[ C os20 2 _ - cos20!_] + -A- [ cos 2 (0 2+ - 2 -) ~ cos2(0 1+ - 0x_)],(A38) 



2AK 



ni y - n 2y = -^-[sin202_ - sin 20^] + — [sin2(0 2+ - 2 _) - sin2(0 1+ - 0i_)](A39) 



We can simplify these expressions by first noting that 

9 a+ - a _ ~ Kt a (l + 0(A/e ± ) 2 ) , (a = 1, 2) 



(A40) 
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and then denning sum and difference variables 



t — \{t\ + t 2 ), tl2 — t\ — £2 



(A41) 



in terms of which we have identities like 



cos 2^2- — cos 26*i_ = 2 sinft_tsin ift_t 



12- 



sin 26*2- — sin 2#i_ 



-2 cos ft_t sin 



12- 



(A42) 
(A43) 



etc. Putting all these together, we obtain 

AAK 



2A 



ni x — n 2x = sinft_t sin ^Q_t 12 H sin 2Kt sin Kt\ 2y 



4AK — 1 2A 

— — cos ft_t sin 7rfi-ti2 — 7^— cos 2Kt sin Kt 
il,il_ z ill 



12- 



(A44) 
(A45) 



Squaring and adding these two expressions, we obtain the mismatch as 



7] = 2 



' AK 



.2 i( 

2 V 



sin" -^-£12 



1 (A 

2 I ft. 



sin 2 



12 



A 2 K - . 

2 —7; - cosft+t sin ift_£ 12 sinKti 2 . 



ftlft_ 



12- 

(A46) 



We can also write this in a manifestly positive form: 



A* 

2ft 2 



K ft_t 12 

2 sin sin Kt\ 2 

ft_ 2 



A 2 A" 

+ 4 sin 2 \VL + t sin \VL_t 12 sin AT£ 12 . (A47) 



Since this expression is unchanged when the signs of both t\ and t 2 are reversed, we have 
now proven our claim that it is also the mismatch for Equations (IA460 and (IA470 are 
valid for t\ ]2 <C E dm /A 2 on account of Eq. (1A40|) . 

To find the phase <j) i n Eq. (IA4I) . we note that from Eqs. (IA29I) and f)A40|) that 



A21 = e-* Kt * + °^ 2 



0(Af 



and likewise for A^. Now, since A a i = 0(A), 



(fii|fi 2 ) = A*^A 21 + A* n A 2i 
= ^ T A 2T + 0(A) 2 
= e lKt12 + 0(A) 2 . 



(A48) 



(A49) 
(A50) 
(A51) 



Thus, the dominant term in is just the zeroth order one, i.e., ~ Kt\ 2 , and we have, 



(A52) 
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with f] given by Eq. (1A46j) . As a check, note that correctly vanishes when tu = 0. 

The calculation above assumes that K ~ e ^> A. We shall see in Appendix [B] that distant 
spins for which K ~ A play an important role in determining the net influence factor. It is 
therefore desirable to find rj when K is small. This can be done by evaluating Eq. (13.91) for 
F^ by a standard perturbation expansion in K. The result is 

K 2 e 2 A 2 K 2 e 2 A 2 K 2 A 2 K 2 A 4 

V = 2Q 4 ^ 2 Q^ tl2 Sin Qtu + 2 ~Q r ~ Sin2 ^ tl2 ~ ~2£$~ Sin2 Qtl2 ' ^ A53 ^ 

Here, f2 = (e 2 + A 2 ) 1 / 2 . The last term is smaller than the first three by order (A/Ed m ) 2 . The 
remaining three terms are qualitatively very similar to what we get from Eq. (1A46I) when 
K -f 0. 



APPENDIX B: ESTIMATE OF MULTI-SPIN INFLUENCE FACTOR FOR 
MOLECULAR SPIN ENVIRONMENT 

In this appendix, we will estimate the total influence factor F = Hi F% f° r I ^12 1 ^ Ej^, 
using a mix of analytic and numerical approaches. 

1. Preliminary Analytic Estimate 

Since 77^ <C 1, the total influence factor is given by 

p ~ e'Ei^e-Ei*, (Bl) 

so \F\ ~ exp(— J2iVi)- We therefore focus on the sum Y^iVi- Let us divide it into three 
parts, Si, 5*2, and 5*3, corresponding to the three terms in r\i in Eq. (|A46|) . 
The first sum is 

For large | i 12 1 (but with |^X2 1 A -1 ), we may replace the factor Q~ 2 sin 2 (f2j_ti2/2) by 
a term proportional to 5(Oj_)|ii2| as in textbook derivations of Fermi's golden rule. The 
replacement must be done with care, however. The physical point is that for large \tu\, 
the only sites that contribute significantly to Si are those for which Qj_ is very small. By 
taking the distribution as 5(fij_), we get a vanishing answer for Si since > A, and so 
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the argument of the ^-function is never satisfied. The correct result which preserves the 
integral with respect to fli- is 



sin(fij_ti2/2) 

a - 



1 2 



^0V-A)|* 12 |. (B3) 

In this equation, the weight of the delta-function is 7r/4 instead of 7r/2 since on the left we 
only integrate over positive values of f2j_, but on the right we wish to interpret the delta- 
function in the standard way, that distribution to be integrated over all f2j_. Using 
Eq. ([B3]) in Eq. flB2J) yields 

Si = l E A2 ^(^-- A )l*i2i- (B4) 

Z i 1 L i+ 

To further simplify this result, we note that = A implies ej_ = 0, e i+ = 2Ki, and 
tt 2 i+ = A 2 + AKf. Therefore, 

5(Oi_ - A) = ^(e* - Ki) = Atffc - Ki), (B5) 

and 

^ = ^A 3 |t 12 | £ J 2 K ^ 2 5{e t - Jti). (B6) 

We now average over the bias distribution (13.41) . This turns 5(ei — Ki) into jiKi)- It then 
remains to do the sum over the sites. Because the summand is slowly varying, we may 
replace the sum by an integral. This integral may in turn be performed by introducing the 
density of couplings g(K), defined so that g(K)dK is the number of sites for which Ki lies 
between K and K + dK. In this way we get 

Sl = r 4 ^ Lc. dK irfe f{K )g{K >■ (B7) 

We have cut off the K integration so as to exclude very distant spins for which the coupling 
is weaker than cA, where c is some constant of order unity. The reason is that for such 
spins the mismatch will be essentially zero, since they are insensitive to the orientation of 
the central spin. 

We show in Appendix [C] that 

g(K) = (B8) 
Note that couplings +K and — K are equally likely. Using this result, we obtain 

Sl = TV T~5T ' tl2 ' L dK K(4K* + A 2 ) 6 K (B9) 
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The integral is dominated by small values of K close to A, so it may be evaluated by parts. 
Doing so, and setting e~ A2 ^ 2E b ~ 1, we obtain 



4tt /2tt . /. 1 \ AE, 



Note that the Gaussian form of /(e) is not essential to the form of the answer. 
The second term in Eq. (1A46j) leads to the sum 

^ = ^E^sin 2 (i^ 12 ). (Bll) 

Z i 1 L i+ 

Since sin 2 {Kitvz) is bounded by 1 the dominant contribution will come from sites on which 
w 0(A 2 ). Averaging over the bias field distribution gives 

1 roc A 2 

52 = 2 ? L ^TJ7TW m siD(K ' ta) *' (B12) 

Because A <C E^m, the integral is very sharply peaked at e = — -fQ. We may therefore 
replace /(e) by the constant /(— Kj). The integral is then elementary, and we obtain 

S 2 = i7rA^/(-^)sin 2 (^ 12 ). (B13) 

i 

The sum is now evaluated as before, by converting to an integral over K. Using the result 
flBEJ) for g{K), we obtain 

This time the integrand is sufficiently convergent near K = 0, so the limit cA can be replaced 
by 0. The factor K~ 2 sin 2 (Ktu) behaves like n\ti2\5(K) for large |^aa | - We may therefore 
replace f(—K) by /(0) = (2ttE 2 )~ 1 ^ 2 , after which the integral is trivial, and yields 

4tt 2 [2^AE dm 

S2 = — VT^r |ti2 '- (B15) 

This result is also valid only for |£ 12 | <C A -1 . 

The last sum, S3, from the third term in Eq. (IA46I) . is given by 



A 2 K 

S z = 2 E TyT-FT sin(fii_ti 2 /2) sin(iQt 12 ) cos(f2 i+ t). (B16) 

As |ti 2 | increases, the term fi^ 1 sin(f2j_t 12 /2) behaves like a 5-function of f2j_. By the same 
reasoning as for S\, we find that the correct replacement is 

sin(£VW 2 ) 7T 

^ = - A ) sgn(ti 2 ). (B17) 
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Further writing Qi± in terms of A, e*, and Ki, we obtain 



3 



& = 4 ^ 2 A + A2 8in(|^t 12 |) cos(^4K 2 + A 2 i) 5(e t - K t ), (B18) 

where we have incorporated the sgn(ti 2 ) and sgni^ factors by taking an absolute value of 
the argument in sin(|i^jt 12 |). The next step is to average over the bias distribution and 
integrate over the sites. As in the case of Si, we exclude distant spins, and obtain, 



S ° = W^^t L iK K^ + A 2 ) *W»I> (B19 ) 
This integral is also dominated by the lower limit, but the answer is different depending on i. 
If t <C A , we may argue that for K ~ A, and for < |i 12 | < A^ 1 , sm(K\t 12 \) « #1*12 1, 
and cos(\/4-ft' 2 + A 2 t) « 1. The resulting integral is identical to that which appeared in Si. 
Hence, we have 

S, = *lJ*u( 1 +l)**l M , (f«A-'). (B20) 



9 V 3 V 4c 2 / E> 



'A, 



If on the other hand t > A x , the oscillations in the cos(\fAK 2 + A 2 i) factor reduce S3 
significantly. The precise from is unimportant, and it suffices to put 

S 3 ~0, (t^A- 1 ). (B21) 

For short times i <C A -1 (which automatically implies \tiz\ <C A -1 ), all three sums Si, 
have the same behaviour. Adding them together, we obtain, 



7r_1M L + 4^ 



A£<,V m . . . 

ti 2 | (B22) 



On the other hand, for longer t, but still obeying A|ti 2 | <C 1, S3 may be neglected, and 

4/T /27T 



9 V 3 



(B23) 



2. Improved estimate incorporating numerics 

Since the analytical estimate given above entails several approximations, we have also 
evaluated \F\ numerically. We take the contribution F,i from the zth MS to be given by 
e lSiKtl2 {\ — with rji given by Eq. (1A46I) . The different factors 77, are found, and the 
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factors (1 — r]i) multiplied to obtain \F\. This calculation is valid for a much larger range of 
times, ti t 2 ~ Ed m /A 2 , since the expression (IA46I) then holds. 

In more detail, our algorithm is as follows. We first create a set of sites on a nearly 
cubic lattice. That is, each site is offset from a perfect cubic lattice by a small random 
amount equal to 0.01 times the lattice constant in each of the three cartesian directions. 
(The reason for adding the offsets was to avoid exact cancellations of the dipole field from 
an aligned shell of nearest neighbor spins. We do not believe that this step is essential, but 
it does not invalidate the calculation either.) We next place a spin on each site with the 
orientation si randomly chosen to be ±1 with equal probability, and choose a particular 
value of ti2- Next, at each site we select an energy bias by sampling a normal distribution 
with mean zero and standard deviation = 1000 in units such that A = 1 . (In Fe§ the ratio 
Eb/A is ~ 10 6 . Taking such a large ratio in the numerics makes each individual mismatch 
prohibitively small, making it very hard to see departures from unity in Y\(l — r]i). The 
physically important point is to ensure Eb/A 3> 1, which we do.) The dipole field Ki at 
each site due to the central spin is computed using Eq. CI3 . 31) with Ed m also equal to 1000A. 
That is, we do not take Ejj and Edm to be different. With these values of Ki and e«, we can 
then find (all energies are computed in units of A) 

e i± = 6i± Q i± 2 = e i± 2 + A 2 . (B24) 

It is now possible to calculate the expression (IA46P for r\i for any t and t^. The dependence 
on two time variables is inconvenient, however, so instead we calculate the lower and upper 
bounds with respect to t, which are given by 
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2 r 



Vi max 2 



A\ 2 



Since |F| = Hi |1 — Vi\i we have 



Q-) *n[-2-)-z*n(Kit*j 

K \ . f Qj-t 12 \ , 1 . ftr . , 
sm — - — + - sin(Ajti 2 J 



Qi-I V 2 J 2 



(B25) 
(B26) 



LKi-^max) < W\ ^n^-^min). (B27) 

i i 

These bounds, |-F| m i n and |-F| max , are now found using the computed maxima and minima 
for rji. At the same time, we also compute the sums Si, 5*2, and 

S 'z = 2E?^7r isin (^i2/2)sin(K i t 12 ), (B28) 



n 2 q 
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which differs from S3 in that the factor cos(f2j+t) is lacking in the summand. As argued 
above, S' 3 ~ S3 when At <C 1. In this time range, therefore, we expect \F\ ~ |-F| m ax- For 
longer i on the other hand, \F\ ~ exp(— (Si + S 2 )). We nevertheless shall find it useful to 
continue to calculate S' 3 . 

We then recompute the S"s, and the bounds for \F\ using a different set of biases, All 
told, we do this for about 10 5 bias configurations, in order to generate averages for the S"s, 
and |F| min and |-F| max . We also calculate the variances in these quantities at this stage. The 
entire calculation is then repeated for different £ 12 . 

The lower and upper bounds of \F\ , |-F| min and |F| max , are plotted as a function of ti 2 
in Figs. US (We do not show the data for N = 8000 and A/E dm values of 0.003 or 0.005.) 

We see that |-F| mm does die as an exponential, i.e., we can fit it to a form e _a l A * 12 l very 
well. On the other hand, |-F| ma x dies like e _fe ' A * 12 l 2 , which is rather different. We show the 
best fit values of a and b for three different A/Edm m Table [H As can be seen, a and b are 
reasonably independent of this ratio. 

The result |-F| max ~ e -b( At i2) 2 j s ra ther surprising, since 



-(s 1+ s 2+ sy (B29) 
e - ( s 1+ s 2 -sy (B30) 



and we showed that all three sums vary linearly with | i 12 1 - We can understand the t\ 2 
behaviour from our numerics. We first note that we do indeed find an excellent linear |ti2| 
variation for the the individual S"s. We therefore write Si = QA\ti2, and determine Q from 
our data. These values are shown in Table HT1 Examining the table, we see that (a) (1 ~ £2? 
and (b) there is a nearly total cancellation in Ci + C2 — (3- i- e -, C3 — Ci + C2- 

We can use the numerical results to improve our analytical estimate as follows. We assume 
that the above mentioned cancellation is perfect. In other words, the unknown parameter c 
in Eq. flB22l is such that 



+ ^) r (1331) 



This implies that 

£77™ = Si + S 2 + S' 3 = —J— [ 
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16vr 2 \2nAE, 
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^^\tn\ (B32) 
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' hn 'ti 2 \. (B33) 



With Edm = Eb, this equals 25.4Ati 2 . Our numerical fits to |F|mm yield a ~ 24, which is 
quite close. This gives us confidence in fixing c as per Eq. (IB31I) [32]. 

We can now estimate J2iVi f° r l° n g t (but <C Ej m /A). Using Eq. (lB3ip . we have 



8tt 2 MkAE, 



dm 



12 



9 V 3 E b 

Since we do not know Ed m /Eb precisely, however, we limit ourselves to stating that 



J^Vi — lrrA\t 



12 



(B34) 



(B35) 



where 7 m is a constant of order unity. 



APPENDIX C: DENSITY OF DIPOLE COUPLING STRENGTHS 

The density of dipole couplings, g(K), introduced in Eq. (1B8|) . is given by 

g(K)=J2S(K-Ki), (CI) 

i 

with 

^ = ^%^(l-3cos 2 ^). (C2) 

' i 

We evaluate the sum over lattice sites assuming that the spins are uniformly distributed with 
a density a~ 3 . Except when K ~ Ed m , corresponding to nearest or next-nearest neighbour 
sites, we may replace the sum by an integral, obtaining 



2yj- roc rl 

g(K) — — - / drr 2 \ du< 
a 6 Jo J-i 



iT-^§^(l-3 M 2 ) 



(C3) 



where u = cos 8. Performing the r integral, we get 



9W = ^/ 1 d „| 1 _ 3 ^|efi^V (04) 



3K 2 7-i 1 1 \ K 

where O(-) is the Heavyside step function; equal to 1 when its argument is positive, and 
zero otherwise. The integral on u is best done separately for positive and negative K. When 
K > 0, we have 

9(K) = ^rJ d«(l-3«) = ^-^-. (C5) 
Likewise, when K < 0, we have 

(V\ 87lE d™ f 1 A (Q 2 -,\ 167r Edm , na , 
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the same expression as for K > 0. This is Eq. ( 1B8I) . 
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When we make these changes, and fix c according to Eq. (|B31|) . we find S1 + S2 — S3 = 6| A.ti2 1 2 , 
with b = 2.02E ( [ m /Ei ) , or just 2.02 if E^m = Ef,. Our numerical fits yield b ~ 2.1-2.8, which is 
close enough given the crudeness of our estimates. 
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TABLE I: Best fit values of the parameters a and b. 







a 






b 






N = 8000 




N = 27000 


N = 8000 




N = 27000 


0.001 


23.7 




24.4 


2.09 




2.80 


0.003 


24.0 




24.2 


2.52 




2.65 


0.005 


23.1 




23.3 


2.64 




2.77 


TABLE II: Numerically calculated values of the coefficients Ci 


, C2, and C3 in the sums 


Si, S*2, and 


S's- 


















N = 8000 






N = 27000 






Ci 


C2 




Ci 


C2 




0.001 


5.44 


5.41 


10.42 


6.06 


6.06 


11.82 


0.003 


6.08 


6.08 


11.90 


6.12 


6.11 


11.77 


0.005 


5.98 


5.99 


11.78 


6.03 


6.11 


11.01 
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FIG. 1: Inferred values of tunnel splitting as a function of the rate at which the applied field is 
swept, assuming that the spin flip probability is given by the Landau-Zener-Stuckelberg formula, 
Eq. (|5.7p . The curves marked K and MA are obtained when the true flip probability is taken 
to obey Kayanuma's formula, Eq. (|5.6p . and the macroscopically averaged formula obtained by 
integrating Eq. (|5.17 ) 



with Fig. 7 of Ref. [5]. 



For the latter, we took 4aEd m /W = 40. This figure should be compared 
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At 12 

FIG. 2: Numerically computed lower and upper bounds |-F| m in and |-F| max , plotted vs. At±2, for a 

central spin in a lattice of 8000 spins. We have chosen A = 0.001-Ed m . The curves are best fits to 
e -a|At la | for | F | m . n) and e -6|At 12 p for | F | max> 
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